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FOREWORD 


The  model  development  work  reported  herein  is  in  support  of  Task 
11860340006  of  the  Geographic  Sciences  Division,  U.  S.  Array  Engineer 
Topographic  laboratories  (ETL),  project  entitled  "Military  Geographic 
intelligence  (MGl)  Products  to  Support  Battlefield  Sensor  Activities." 

The  objective  o*”  the  project  is  to  design  a  family  of  prototype  military 
geographic  intelligence  products  to  support  planning  for  the  use  and 
operational  placement  of  ground-contact  sensors  on  a  battlefield  to 
detect  the  presence  of  enemy  troops  and  equipment. 

The  U.  S.  Army  Engineer  Waterways  Experiment  Station  (wSS)  con¬ 
tribution  to  the  MGI  project  depended  heavily  on  data  collected  in  a 
number  of  seismic  sensor  programs.  Acknowledgment  is  given  for  data 
furnished  by  Project  MASSTER,  Defense  Special  Projects  Group  (DSf.FG), 
the  U.  S.  Army  Test  and  Evaluation  Command,  and  the  Office,  Chief  of 
Engineers . 

The  work  reported  herein  is  the  result  of  a  coordinated  effort 
during  the  period  1  January-30  June  1971  by  members  of  the  Terrain 
Analysis  Branch,  Mobility  and  Environmental  (M&E)  Division;  the  Soil 
Dynamics  and  Geology  Branches,  Soils  Division;  and  the  Operations 
Branch,  Instrumentation  Services  Division.  Key  participants  in  the 
study  were  Messrs.  Bob  0.  Benn  and  L.  E.  Link  of  the  Terrain  Analysis 
Branch  and  Mr.  Robert  F.  Ballard  and  Dr.  William  F.  Marcuson  of  the 
Soil  Dynamics  Branch.  The  report  was  prepared  by  Messrs.  Benn  and  Link. 

The  study  was  under  the  direct  supervision  of  Mr.  Benn,  the  Pro¬ 
gram  Manager,  and  under  the  general  supervision  of  Mr.  W.  E.  Grabau, 

Chief,  TAB,  and  Messrs.  W.  G.  Shockley  and  S.  J.  Knight,  Chief  and 
Assistant  Chief,  respectively,  of  the  M&E  Division;  The  Director  of 
the  WES  during  the  study  period  was  COL  Ernest  D.  Peixotto,  CE.  The 
Technical  Director  was  Mr.  F.  R.  Brown. 
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SUMMARY 


Improved  guidance  manuals  Tor  planning  the  deployment  and  emplace¬ 
ment  of  seismic  intrusion  detectors  (SID's)  are  needed  to  optimize  the 
use  of  these  devices  for  battlefield  surveillance-  The  development  of 
these  Military  Geographic  Intelligence  (MGl)  products  requires  a  de¬ 
tailed  understanding  of  the  operating  principles  of  the  detector  coupled 
with  an  equally  detailed  understanding  of  the  interactions  of  the  sensor 
propagation  mode  with  the  operational  environment.  This  report  presents 
the  results  of  a  preliminary  analysis  of  data  collected  in  a  wide  range 
of  environments  at  22  sites  in  Panama,  10  sites  in  "Puerto  Rico,  6  sites 
near  Yuma  Proving  Ground,  Arizona,  and  9  sites  near  Ft.  Huachuea, 
Arizona. - 


Multiple  regression  techniques  were  used  to  determine  the  terrain 
factors  that  could  be  correlated  with  the  seismic  responses  resulting 
from  a  man  walking  or  a  controlled  source  (drop  hammer)  that  simulated 
the  signature  resulting  from  a  footstep^  The  measure  of  seismic  re¬ 
sponse  was  peak  particle  velocity  as  a  function  of  distance  from  the 
source.  The. -terrain  factors  that  correlated  best  with  peak  particle 
velocity  were  the  thickness  of  the  first  refraction  layer,  cone  index 
of  the  Or  to  15-*cm  soil  layer,  dry  density  of  surface  soil  and  first  soil 
layer,  water  content  of  surface  soil  and  first  soil  layer,  compression 
wave  Velocity,  Rayleigh  wave  velocity,  and  grain-size  distribution.  An 
empirical  equation  was .developed  to  predict  peak;  particle  velocity  ver¬ 
sus  distance  as  a  function  of  the  terrain,  factors . 

Hies  required  to  trigger  the  logic  of  the  Phase 
imposed  on  the  predicted  peak  particle  velocity 
prediction  of  sensor  performance.  These  computation 
computerized  to  make  a  prediction  model  for  relative 
a  function  of  terrain  factor  values. 

The  empirical  prediction  equation  adequately  predicted 
particle  velocity-distance  relation;  however,  the  predictions 
performance  were  inadequate .  The  errors  in  the  predictions  of 
performance  were  attributed  to  the  inadequacy  of  the  peak  particle 
velocity-distance  relation  to  represent  the  complex  interaction  of  the 
entire  seismic  signal  and  the  sensor.  Frequency  characteristics  of  the 
seismic  signal  and  the  frequency  response  characteristics  of  the  sensors 
also  must  be  considered.  .  ' 
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EFFECTS  OF  ENVIRONMENT  ON  SEISMIC  INTRUSION 
DETECT®  PERFORMANCE 


PART  I:  INTRODUCTION 
Background 

1.  The  urgent  need  for  battlefield  surveillance  devices  has 
prompted  the  rapid  development  bf  seismic  and  other  sensor  hardware  that 
have  proven  useful  in  a  number  of  combat  situations  in  South  Vietnam. 
Unfortunately,  success  stories  cannot  be  written  for  all  sensor  deploy¬ 
ment  attempts.  A  major  factor  recognized  as  responsible  for  the  less 
than  outstanding  performance  of  the  devices  (Military  Geographic  Intel¬ 
ligence  (MGl)  products)  is  the  lack  of  adequate  guidance  manuals  for 
planning  their  deployment  and  emplacement.  Attempts  to  fill  this  gap 
have  not  been  successful  because  the  production  of  such  manuals  requires 
a  detailed  under standing  of  the  operating  principles  of  the  detectors 
coupled  with  an  equally  detailed  understanding  of  the  interactions  of 
the  sensor  propagation  mode  with  the  operational  environment.  Consid¬ 
erable  development  work  is  known  to  be  needed  before  rational  procedures 
for  producing  MGI  products  for  all  types  of  sensor  systems  can  be  formu¬ 
lated.  For  example ,  although  this  effort  has  been  restricted  to  seismic 
intrusion  detectors  (SID’s,),  the  results  have  to  be  considered  interim 
solutions  until  more  definitive  work  can  be  completed. 

2.  The  operating  principles-,  of  SID’s  are  reasonably  well  known; 
however,  the  manner  in  which  the  seismic  energy  is  transferred  from  the 
source  (vehicle  or  persons)  to  the  ground,  the  way  the  substrate  con¬ 
ditions  affect  the  energy  propagation,  and  the  manner  in  which  the  en¬ 
ergy  is  transferred  from  the  ground  to  the  sensor  geophone  are  not  un¬ 
derstood.  Theoretical  solutions  to  these  problems  are  and  have  been  the 
partial  objectives  of  considerable  research  sponsored  by  the  Department 
of  Defense.  The  results  of  current  research  in  this  a’ea  being  con¬ 
ducted  at  the  U.  S.  Array  Engineer  Waterways  Experiment  Station  (WES)  for 
the  Office,  Chief  of  Engineers  (OCE),  and  the  Defense  Special  Projects 
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Group  (DSFG)  are  encouraging;  nevertheless,  a  practical  theoretical 
solution  will  require  additional  research. 

3.  Work  on  the  projects  mentioned  above  and  other  related  efforts 
have  resulted  in  "the  collection  of  considerable  data  that  are  available 
for  empirical  analysis.  Some  analysis  has  been  accomplished  and  an  equa¬ 
tion  has  been  formed  that  allows  prediction  of  seismic  signal  levels  as  a 
function  of  distance  from  a  source.  This  equation  has  been  coupled  with 
seismic  sensor  performance  specifications  to  formulate  an  interim  sensor 
performance  prediction  model. 

4.  Specific  and  quantitatively  defined  terrain  factors  are  the 
inputs  to  the  interim  model.  To  obtain  these  inputs  easily  from  conven¬ 
tional  terrain  intelligence  gathering  techniques  would  be  desirable;  but 
not  all  the  present  terrain  inputs  can  be  so  obtained,  and  transforms 
must  be  found. 

5 .,  .  The  equation  used,  in  the  interim  prediction  model  has  been 
derived  from  the  analysis  of  data  collected  in  a  wide  range  of  environ¬ 
ments,  i.e.  Puerto  Rico,  Panama,  and  Arizona.  Additional  data  will  be 
collected  and  the  prediction  model  improved  as  the  data  are  analyzed. 

The  interim  model  has  been  designed  so  that  modifications  can  be  made 
easily  as  new  data-  or  theoretical  information  becomes  available. 


6.  The  purpose  of  the  study  reported  herein  was  to  provide  techni¬ 
cal  assistance  to  the  U.  S.  Army  Engineer  Topographic  Laboratories  (-ETL) 
seismic  sensor  MGI  product  project  by  developing  theoretical  relations 

f 

of  seismic  wave  propagation  in  .earth  media  and  empirical  relations  based 
on.  the  analysis  of  existing  terrain/sehsor  data. 

7.  The;  report,  includes,  a  brief  discussion  of  the  seismic  and  en¬ 
vironmental  field  -data  collection  programs,,  and  data  reduction  and  anal¬ 
ysis  procedures 'followed  in  the  development  pf  a  seismic  sensor  perfor¬ 
mance  model.  Techniques  for  measuring  or 'estimating  the  terrain  factor 
inputs  to  the  model  are  also  discussed. 

8.  Existing  theoretical  solutions  were  found  to  be  in  a  form, not 
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directly  compatible  with  the  project  objective.  They  were  used  exten¬ 
sively,  however^  in  the  design  of  the  data  collection,  reduction,  and 
analysis  procedures  used  in  the  development  of  the  model. 

Definitions 

9.  Certain  terms  pertinent  to  this  study  and  having  restricted 
meaning  are  defined  below. 

Cone  index.  ^  An  index  of  the  shearing  resistance  of  a  medium  ob¬ 
tained  with  the  cone  penetrometer.  The  value  represents  the  resistance 
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of  the  medium  to  penetration  of  a  30-deg  cone  of  0.5-sq-in.  (6.45-em  ) 
base  or  projected  area.  The  number,  although  usually  considered  dimen¬ 
sionless  in  trafficability  studies,  actually  denotes  pounds  of  force  on 
the  handle  divided  by  the  area  of  the  cone  base  in  square  inches.  The 
cone  index  of  the  soil  surface  and'  the  average  cone  index  of  the  0-  to 
15- cm  layer  are  used  in  this  report. 

Dry  density  (74).^  Dry  unit  weight;  the  weight  of  oven-dried  soil 
solids  (Ws )  from  a  sample  per  unit  of  total  volume  (VT)  of  the  soil  sam¬ 
ple.  Symbolically  this  is 

Ws  - 

7  =-=■( in  g/cnr  ) 

T 

Particle  velocity.  The  time  rate-  of  change  of  the  motion  of  a 

particle  of  the  medium  with  respect  to  a  specified  reference  frame.  The 

particle  velocity  was  measured  by  the  geophones  used  in  this  study. 

Water' content  (w)  .1  The  ratio  of  the  weight  of  water  (W  ')  in  a 

sample  of  soil  to  the  weight  of  soil  (solids  only)  (W  )  in  the  same  sam- 

s 

pie  expressed  as  a  decimal.  It  may  be  written  as 


Compression  wave  velocity  (-Vp).-  The  speed  of  a  compression^  wave 
through  a  medium.  Compression  waves  have  the  greatest  velocity  of  any 
elastic  wave  in  the  same  medium.-  The  motion  of  the  particles  is 


parallel  to  the  direction  of  propagation.  Vp  is  defined  mathematical 1 y 
as- 


where 

V-,  =  compression  wave'  velocity,  LT  ■*" 

F  -1-2 
\  =  Lame's  constant,  ML  T 

-2 

G  =  shear  modulus,  FL 

-4  2 

p  =  mass  density,  GL  T 

Rayleigh  wave  velocity  (Vr).  The  speed  of  a  Rayleigh  wave  (par¬ 
ticle  motion  is  elliptically  retrograde  and  parallel  to  the  direction  of 
propagation)  along  the  free  surface  of  a  medium;  depends  on  Poisson's 
ratio  (v)  of  the  medium.  For  values  of  Poisson's  ratio  0<v<0.5,  the 

Rayleigh  wave  velocity  has  the  range  0.875V-  <v  <0.955V_,  where  V„ 

b  T  b  b 

=  shear  wave  velocity. 

Shear  wave  velocity  (Vg).  The  speed  of  a  shear  wave  (particle 
motion  of  the  medium  is  perpendicular  to  the  direction  of  propagation) 
thrpugh  a  medium,  defined  mathematically  by  the  equation 


where 

-1’ 

V  =  shear  wave  velocity,  LT 

b  -2 

G  =  shear  modulus,  FL 

.  -4  2 

p  =  mass  density,  GL  T 

Thickness'  of  the  first  soil  layer  (%).  The  vertical  depth  (i.e. 
perpendicular  to  the  surface)  to  the  interface  between  the  surface  layer 
and  the  next  shallowest  .layer  as  distinguished  by  their  differing  pri¬ 
mary  wave  velocities.  The  primary  wave  velocities  of  these  two  layers 
are  determined  by  techniques  of  refraction  seismology.  (Note:  The 
above-defined  layers  often,  but  not  necessarily,  correspond  to  soil 
layers  as  defined  ‘by  nonseismic  parameters  (e.g.  grain  size,,  density. 


PART  II:  DERIVATION  OF  TERRAIN/SENSOR  INTERACTIONS 


Field  Data  Collection  Program 

10.  The  ETL  seismic  sensor  study  utilized  data  collected  in  the  con¬ 
duct  of  related  seismic  sensor  research.  The  field  data  collection  pro¬ 
gram  was  designed  to  provide  information  for  development  of  empirical 
terraih/seismic  response  relations  and  to  verify  theoretically  developed 
relations.  The  approach  used  was  to  perform  special  seismic  tests  in 
various  environmental  conditions,  collect  environmental  data  concur¬ 
rently  with,  the  seismic  tests,  and  then  study  these  data  to  determine 
empirically  the  effects  of  environment  on  seismic  response.  This  part 

of  the  report  discusses  the  field  data  collection  program,  derivation  of 

particle  velocity/environment/distance  relations  and  the  maimer  in  which 
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the  relations  were  coupled  with  Phase  III  sensor  logic  to  provide  a 
capability  for  predicting  an  indicator  of  sensor  performance.  Phase 
III  sensors  are  the  most  recently  developed  SID's;  the  sensor  develop¬ 
ment  program  be -an  with  the  Phase  I  sensors  and  has  since  moved  through 
the  Phase  II  sensors  into  the  development  of  improved  devices  that  are 
designated  Phase  III  sensors. 

Site  selection  ; 

11.  To  ensure  that  tests  were  conducted  in  a  wide  range  of  en¬ 
vironmental  conditions ,  care  was  taken  in  selecting  sites  within  pre¬ 
determined  study  areas.  Field  work  was  accomplished  at  22  sites  in 
Panama,  10  sites-  in  Puerto  Rico,  and  15  sites  in  Arizona.  The  sites 
were  tentatively  selected  in  the  office  by  utilization  of  available 
published  data,  topographic  maps,  and  air  photos.  The  published  data 
and  maps  were  used  to  supplement  a  photo  interpretation  study  that 
involved  a  stereoscopic  examination  of  the  photos  of  the  study  area. 

Photo  patterns  were  isolated  on  the  basis  of  their  tone,  texture,  and 
shape;  and  the  assumption  was  made  that  each  discrete  pattern  repre¬ 
sented  a  certain  combination  of  environmental  conditions.  Sites  were 
tentatively  selected  to  encompass  as  many  of  the  terrain  conditions  as 
possible  within  a  study  area;  however,-  accessibility  to  a  site 
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was  also  taken  into  consideration. 

12.  The  final  selection  of  sites  was  not  made  until  after  a 
ground  reconnaissance  was  performed  to  determine  the  validity  of  the 
assumptions  made  during  the  office  study.  Soil,  surface  geometry,  and 
vegetation  conditions  were,  observed  at  the  tentative  sites.  Other 
areas  were  also  visited  to  see  if  different  conditions  existed  that  were 
not  recognized  during  the  photo  interpretation  study.  Upon  completion 
of  the  ground  reconnaissance,  the  final  site  selections  were  made. 

Site  layout 

13.  To  perfdrm  the  seismic  tests,  each  site  was  prepared  in  a 
specific  manner.  A  walk  path  was  laid  out  60  m  long,  with  the  0-  sta¬ 
tion  located  in  the  middle  and  stakes  at  5-m  intervals  along  the 
entire  length  of  the  path.  Another  line  of  the  same  length  was  laid 
out  perpendicular  to  the  first  and  intersecting  it  2  m  from  the  0- 
station.  Stakes  also  were  placed  at  5-m  intervals  along  the  second  line. 
Fig.  .1  shows  a  typical  site  layout. 

Environmental  data 

14.  All  of  the  environmental  factors  that  were  hypothesized  to 
affect  the  seismic  response-  of  the  area  were  considered.  The  fol¬ 
lowing  environmental  descriptors  were  believed  to  be  important,  and 
information  was  collected  on  each  either  by  direct  measurement,  labora¬ 
tory  analysis,  or  computation. 

a.  Soil  characteristics 

(1)  Layer  thickness  (refraction  seismic  technique) 

(2)  Moisture  -content 

(3)  Dry  density 
(4:)  Void  ratio 

c  (5)  Degree  of' saturation 

(6)  Liquid  limit 

(7)  Plastic  limit 

(8) -  Cone  index 

(9)  Grain- size  distribution 

b.  Vegetation  characteristics 
(l)  Stem  diameter 


Fig.  1.  Typical  site  layout 


(2)  Stem  height 

(3)  Stem  spacing 

(4)  Crown  volume 

c.  Surface  geometry:  surface  profiles 

d.  Meteorological  conditions 

(1)  Wind  speed 

(2)  Wind  direction 

(3)  Rainfall 

(4)  Air  temperature 

(5)  Soil  temperature 

15-  The  soils  data  were  obtained  by  measuring  the  basic  soil 
parameters  in  the  field  with  such  devices  as  a  nuclear  moisture-density 
meter,  and  by  using  conventional  soil  sampling  procedures  and  obtaining 
values  through  laboratory  analysis  of  the  samples.  The  vegetation  and 
surface  geometry  data  were  collected  according  to  standard  WES  pro¬ 
cedures.  Meteorological  data  were  obtained  with  instrumentation  avail¬ 
able  at  the  test  sites  or  with  a  portable  field  unit  designed  at  the 
WES. 

16.  Techniques  for  measurement  or  calculation  of  the  terrain 

factor  values  have  been  extensively  documented  in  the  references  listed 

at  the  end  of  the  text  of  this  report.  Up-to-date  instrumentation  and 

techniques  were  used,  and  the  data  were  recorded  in  formats  compatible 

with  automatic  data  processing  that  allowed  their  efficient  analysis. 

Complete  documentation  of  the  field  data  collection  program  will  be 
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published  as  part  of  a  -report  dealing  with  a  related  seismic  sensor 
program. 

Seismic  response  data 

17.  The  seismic  response  data  were  collected  at  the  various  loca¬ 
tions  by  measuring  the  particle  velocity  resulting  from  (a)  a  man  walk¬ 
ing  and  (b)  a  controlled  energy  source  (hammer  drop).  The  seismic  re¬ 
sponses  were  measured  with  two,  three-directional  geophones  and  recorded 
on  magnetic  tape  with  a  wide-band  -amplifier-recorder  system.  The 
geophones,  moving-coil  type  with  a -usable  flat  frequency  response  from 
1.5  to  200  Hz,  were  buried  flush  with  the  ground  surface  and  5  m  apart. 
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positioned  as  shown  in  fig.  1.  The  drop-hammer  energy  source  was  de¬ 
signed  with  the  hope  that  its  response  would  be  comparable  in  magnitude 
and  signature  to  the  response  produced  by  a  footstep.  Since  footsteps 
are  very  variable  in  character,  the  particle  velocity  resulting  from 
the  hammer  drop  would  hopefully  allow  a  more  accurate  comparison  of 
seismic  responses  in  various  environmental  conditions. 

18.  In  the  man-walking  tests  a  man  walked  at  a  constant  rate 
along  both  prescribed  paths.  He  started  30  m  from  the  centermost  geo¬ 
phone  array,  and  continued  on  the  same  line  until  he  was  30  m  past  the 
centermost  array,  which  brought  him  to  within  2  m  of  the  centermost 
array  on  each  path.  The  controlled- source  (drop-hammer)  tests  were 
conducted  by  dropping  the  calibrated  weight  of  the  hammer  at  5-m 
intervals  along  the  same  paths  used  for  the  walk  tests. 

19.  The  peak  particle  velocity,  or  maximum  signal  amplitude,  re¬ 
sulting  from  each  footstep  or  hammer  drop  was  obtained  from  the  magnetic 
tape  recordings  by  machine  processing  at  the  WES.  These  data  were  used 
to  develop  peak  particle  velocity-distance  relations.  (The  data  col¬ 
lected  in  Puerto  Rico  were  recorded  on  oscillographs  and  were  reduced 
manually  to  obtain  the  peak  particle  velocity  values  for  the  Puerto 
Rico  sites.) 


Derivation  of  Equations 


Regression  technique 

20.  The  data  collected  in  the  field  in  Panama,  Puerto  Rico,  and 
Arizona  were  used  as  a  base  for  generating  an  empirical  equation  des¬ 
cribing  the  seismic  response.  The  dependent  variable  selected  for  this 
study  was  the  peak  particle  velocity  resulting  from  the  calibrated 
(drop-hammer)  source.  The  independent  variables  consisted  of  the 
various  descriptors  of  the  soil  characteristics  (paragraph  14) ,  com¬ 
binations  of  soil  descriptors,  functions  relating  the  soil  descriptors, 
and  various  combinations  of  all  of  the  above. 

21.  To  formulate  a  peak  particle  velocity-distance  prediction 
capability,  a  basic  format  had  to  be  selected  for  use  in  developing 
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the  empirical  equation.  The  best  results  were  believed  to  be  obtainable 
if  the  empirical  equation  format  conformed  closely  to  theoretical  con¬ 
siderations  of  seismic  energy  decay  with  distance.  Two  types  of  decay 
were  considered,  geometric  damping  and  equivalent  viscous  damping.  Geo¬ 
metric  damping,  the  decay  of  energy  due  to  the  spreading  of  the  wave 
front  over  a  larger  and  larger  volume,  for  a  Rayleigh  wave  can  be 
described  by  1/Vr  ,  where  r  is  the  distance  from  the  source.^ 
Equivalent  viscous  damping  can  be  approximated  by  an  exponential  decay 
function  and  has  been  described  by  the  expression  ^ 


where 

e  =  base  of  natural  logarithm 

k'  =  damping  coefficient 

W  =  mean  circular  frequency 

Vn  =  Rayleigh  wave  velocity 
K 

r  =  radial  distance  from  source 

22.  A  similar  expression  was  chosen  for  the  empirical  analysis: 

Ae~°*r 

Yr 

where  l/Vr  describes  the  geometric  damping  of  the  Rayleigh  wave,  and 
e-Qr  conforms  to  the  equivalent  viscous  damping  of  the  wave,  with 
a  replacing  (k'W)/(2TTVR)  used  in  the  theoretical  expression.  The 
A  term  represents  an  initial  particle  velocity  amplitude  that  is 
attenuated  by  the  e'^/Vr  expression.  In  theoretical  work,  the  A 
term  comprises  a  number  of  functions  that  theoretically  describe  an 
unattenuated  Rayleigh  wave.  In  this  case,  or  is  somewhat  similar  to 
the  equivalent  viscous  damping  coefficient,  although  it  includes  the 
other  terms  in  the  total  coefficient  of  the  theoretical  equation.  For 
the  purpose  of  this  study,  &  was  considered  to  be  a  constant  over  the 
frequency  range  of  the  data  since  a  large  range  in  frequencies  did  not 


occur. 


Description  of  computer  software 

23-  A  multiple  regression  computer  program  written  by  Mr.  J.  H 
Goodnight  of  North  Carolina  State  University,  Raleigh,  North  Carolina, 
was  utilized  to  correlate  the  independent  variables  and  combinations 
thereof  with  the  dependent  variable.  In  general,  use  of  the  computer 
program  entails  four  basic  steps: 

a.  Values  for  the  basic  variables  are  fed  in,  and  values  for 
the  combined  variables  are  generated. 

b.  The  simple  statistics  (i.e.  sum,  mean,  sum  of  squares, 
variance,  and  standard  deviation)  of  each  variable  are 
computed;  and  a  bivariant  analysis  is  conducted,  i.e. 
each  variable  is  correlated  with  every  other  variable, 
one  at  a  time.  This  indicates  the  variables  that  are 
interrelated. 

c.  A  model  is  built  for  each  specified  dependent  variable. 
The  computer  searches  the  independent  variables  (first 
taking  one  at  a  time,  then  two  at  a  time,  and  so  on)  and 
lists  the  individual  variables  that  correlate  best  with 
the  dependent  variable.  These  lists  of  variables  are 
termed  models. 

d.  Based  on  the  models  generated,  the  independent  and  depen¬ 
dent  variables  are  specified,  and  the  computer  uses  a 
Doolittle  matrix  inversion  technique  to  generate  the  re¬ 
gression  equation  of  best  fit  through  the  data.  This 
equation  is  in  the  form 


Y  =  B  +  B.X.  +  B0X0...+  BaXw 
0  xi  2  2  Ifl 

where 

Y  =  dependent  variable 
X^  =  independent  variables 
B^  =  regression  coefficients 

24.  The  correlation  coefficient  R  is  defined  and  used  as  a  mea¬ 
sure  of  "goodness  to  fit."  R  =  1  is  a  perfect  correlation;  whereas 

R  =  0  indicates  no  correlation  at  all.  An  analysis  of  variance  table 
also  is  printed  out  for  determining  the  significance  of  the  equations. 
Use  of  computer  program 

25.  The  relations  between  seismic  response  and  environmental 
factors  have  been  shown  to  be  very  complex.  To  use  the  strength  of  the 
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multiple  regression  technique  to  best  advantage,  a  step-by-step  proce¬ 
dure  for  building  the  desired  relations  from  the  basic  parameters  had  to 
be  followed: 

a.  The  basic  environmental  parameters  and  simple  combinations 
thereof  were  correlated  with  a  given  seismic  response  de¬ 
pendent  variable  to  determine  the  parameters  that  cor¬ 
related  best  with  the  dependent  variable. 

b.  Equations  were  written  for  the  best  models  developed. 

c.  The  equations  were  used  to  determine  whether  the 
variables  were  positive  or  negative  with  respect  to  the 
value  of  the  dependent  variable.  If  the  independent 
variable  appeared  in  a  negative  term  of  the  equation,  it 
was  considered  negative;  if  the  independent  variable 
appeared  in  a  positive  term  in  the  equation,  it  was 
considered  positive. 

d.  The  negative  and  positive  variables  were  then  condensed 
to  generate  new  combined  variables.  Negative  variables 
were  always  in  the  denominator  of  the  new  variables,  or 
as  negative  terms;  whereas  the  positive  variables  were 
always  placed  in  the  numerator  of  the  new  variables, 

or  as  positive  terms. 

e  The  basic  parameters  and  the  newly  generated  combined 
variables  were  again  correlated  with  the  dependent  ■ 
variable  to  improve  the  correlation. 

f.  Equations  were  written  for  the  best  models,  and  the  sta¬ 
tistics  of  the  equations  were  evaluated  to  determine 
whether  further  development  was  required. 

g.  If  further  development  was  required,  a  third  group  of 
combined  variables  were  derived  from  the  newly  written 
equations,  and  the  cycle  was  continued  until  no  addi¬ 
tional  improvements  could  be  made. 

Peak  particle  velocity  relation 

26.  The  relation  between  peak  particle  velocity  and  the  environ¬ 
mental  factors  was  developed  in  two  major  phases,  both  of  which  followed 
the  steps  in  the  procedure  outlined  above.  The  first  phase  was  con¬ 
cerned  with  the  development  of  relations  for  A  and  ot  •  Once  these 
relations  were  determined,  the  total  relation  for  peak  particle  veloc- 
ity,  Ae  /V r  ,  was  developed: 

Ae_ar 

U  peak  =  0.11  +  ~— 

where 
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particle  velocity  resulting  from  the  hammer  drop, 
cn/sec  x  10-3 

radial  distance  of  source  from  the  geophone,  m 
thickness  of  first  soil  layer,  m 
moisture  content  of  the  surface  soil,  percent 
moisture  content  of  the  first  soil  layer,*  percent 

O 

dry  density  of  the  first  soil  layer,*  g/cnr 

q 

dry  density  of  surface  soil,  g/cnr 

cone  index  for  the  0-  to  15-cm  surface  soil  layer,*  psi 
base  of  natural  logarithm 

grains  finer  than  0.074  mm  by  weight  for  first  soil 
layer,*  percent 

compression  wave  velocity  of  first  soil  layer,  n/sec 
Rayleigh  wave  velocity,  n/sec 


27.  A  number  of  standard  statistical  tests  were  used  to  evaluate 
the  empirical  equation.  An  analysis  of  variance  was  performed  for  the 


* 


Average  for  the  layer. 


dependent  variable,  U  peak,  of  the  equation  and  the  resulting  variance 
ratio  F  was  tested  at  the  0.01-  and  0.05-levels  of  significance.  The 
F  test  showed  the  equation  to  be  significant  at  the  0.01  level.  A 
student’s  t-test,  associated  with  the  hypothesis  that  the  regression 
coefficients  are  equal  to- zero,  was  applied  to  the  coefficients  of  the 
regression  equation.  The  results  of  the  t-test  showed  the  coefficients 
to  be  nonzero  and  therefore  statistically  valid  at  the  0.01-significance 
level. 

28.  The  correlation  coefficient-  E  of  the  equation  is  0.83. 

This  indicates  that  the  regression  equation  cam  be  used  to  predict  the 
peak  amplitude  of  a  seismic  signal  with  some  degree  of  precision.  The 

O 

standard  deviation  a  of  the  equation  is  0.77  X  10_J  cn/sec,  which  is 
very  significant  with  respect  to  the  measured  values  of  peak  particle 
velocity  at  distances  greater  than  30  m.  This  indicates  that  the  de¬ 
veloped  equations  will  yield  predictions  that  will  exhibit  considerable 
scatter. 

Prediction  of  Sensor  Performance 


Seismic  sensor  characteristics 

29-  A  schematic  drawing  of  the  major  components,  of  a  seismic 
sensor  is  shown  in  fig.  2.  The  ground  motion  resulting  from  a  seismic 
wave  is  measured  by  the  geophone  and  converted  into  an  analog  electrical 
signal.  The  frequency  and  amplitude  of  the  electrical  signal  are  pro¬ 
portional  to  the  particle  velocity.  The  electrical  signal  from  the 
geophone  is  then  fed  through  a  band-pass  amplifier,  where  the  signal 
is  filtered  and  amplified.  From  the  amplifier  the  signal  goes  to  the 
sensor  logic.  If  the  signal  has  ah  amplitude  above  a  certain  threshold 
value,  the  logic  will  be  activated.  The  logic  will  then  integrate  suc¬ 
ceeding  signals  from  succeeding  footsteps  until  enough  energy  is  com¬ 
piled  to  reach  a  second  threshold,,  which  causes  the  sensor  to  transmit 
a  coded  EF  signal  to  a  receiving  station,  indicating  that  a  source  of 
seismic  signal.'.,  such  as  a  man  walking,  is  nearby. 

30.  It  must  be  emphasized  at  this  point  that  the  signal  reaching 
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Schematic  of  seismic  sensor 


the  logic  is  very  dependent  on  frequency.  The  response  of  the  sensor 
geophone  is  not  constant  for  all  frequencies,  and  the  band-pass  ampli¬ 
fier  attenuates  signals  above  and  below  40  and  0  Hz,  respectively.  Thus 
the  characteristics  of  the  signal  reaching  the  sensor  logic  may  or  may 
not  be  directly  analogous  to  the  frequency  and  amplitude  of  the  ground 
motion. 

31.  For  this  preliminary  study  the  effects  of  signal  frequency 
have  been  ignored  completely,  and  the  signal  amplitude  reaching  the  sen¬ 
sor  logic  has  been  assumed  to  be  equivalent  to  the  peak  particle 
velocity.  The  values  used  for  sensor  logic  thresholds  are  nominal 
values  based  on  design  specifications.  The  actual  threshold  values  in 
the  field  sensors  may  vary  considerably  from  one  sensor  to  another  be¬ 
cause  of  the  wide  tolerance  in  the  manufacturing  specifications  imposed 
to  limit  unit  costs. 

Technique  for  prediction 
of  sensor  performance 

32.  The  empirical  equation  derived  by  the  multiple  regression 
technique  has  been  combined  with  nominal  values  of  sensor  logic  thresh¬ 
olds  to  provide  the  capability  for  predicting  sensor  performance,  as 
she™  in  fig.  3.  The  equation  is  used  to  predict  the  curve  of  peak 


Fig.  3.  Generalized  relations  showing  how  sensor 
performance  can  be  predicted 
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particle  velocity  versus  distance  for  a  given  set  of  environmental  con¬ 
ditions.  Nominal  sensor  logic  threshold  values  for  high,  medium,  and 
low  gain  are  superimposed  on  the  predicted  curve;  the  distance  D  ,  where 
the  threshold  value  for  any  particular  gain  intersects  the  peak  particle 
velocity  curve,  is  the  measure  of  seismic  sensor  performance  for  that 
gain.  This  technique  has  been  computerized  to  provide  an  automated  per¬ 
formance  prediction  capability.  The  computer  program  is  presented  and 
discussed  in  detail  in  Appendix  A. 

Evaluation  of  Prediction  Capability 

33-  The  primary  factor  affecting  the  sensor  performance  predic¬ 
tions  is  the  accuracy  of  the  predicted  particle  velocity-distance  data, 
i.e.  the  adequacy  of  the  regression  equation  discussed  previously. 
Adequacy  of  regression  equation 

34.  The  accuracy  of  the  regression  equation  depends  on  the  close¬ 
ness  of  fit  of  the  regression  equation  to  the  actual  particle  velocity- 
distance  data  measured  in  the  drop-hammer  tests.  Predicted  and  measured 
peak  particle  velocity  data  for  selected  sites  in  Panama  and  Yuma 
Proving  Grounds  are  compared  in  plates  1-3-  In  a  majority  of  cases  the 
predicted  values  compare  closely  with  the  measured  data.  Much  more 
error  occurs  in  the  near-field  portion  of  the  plot  than  in  the  far-field 
portion  (i.e.  at  distances  greater  than  10  m).  These  plots  demonstrate 
that  the  multiple  regression  prediction  equation  is  a  fairly  good  pre¬ 
dictor  of  peak  particle  velocity  for  distances  greater  than  10  m. 
Significance  of  predictions 

35 •  The  computer  program  in  Appendix  A  was  used  to  obtain  pre¬ 
diction  for  medium-  and  low-gain  detection  ranges  for  19  sites  in 
Panama  and  Yuma  Proving  Grounds.  To  evaluate  the  adequacy  of  these 
predictions  for  representing  the  detection  distance  for  a  man  walking, 
a  body  of  reference  data  from  which  comparisons  could  be  made  was 
necessary.  A  portion  of  the  necessary  reference  data  was  available  in 
the  form  of  one  man  walking  in  the  detection  ranges  for  a  DSID  Phase  III 
sensor  at  low  gain  for  each  test  site  in  Panama.  A  total  of  five  tests 
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were  made  at  each  site.  The  data  were  collected  concurrently  with  the 
field  data  collection  program  in  Panama;  however,  they  were  not  con¬ 
sidered  adequate  since  they  contained  a  large  amount  of  scatter  that 
could  not  he  explained.  To  complete  the  necessary  body  of  reference 
data,  detection  distances  were  determined  with  an  analog  computer. 

3 6.  A  model  of  a  Phase  III  seismic  intrusion  detector  with  PID 
logic  was  programmed  on  an  analog  computer  to  form  a  synthetic  sensor, 
and  ■;'he  measured  analog  records  from  the  man-walking  tests  were  used  as 
inputs  to  the  model.  Detection  ranges  for  one  man  walking  were  obtained 
for  both  medium  and  low  gains.  No  attempt  was  made  to  predict  detection 
ranges  associated  with  the  high-gain  setting  because  background  seismic 
noise  levels  at  every  site  thus  far  visited  in  this  program  have  ex¬ 
ceeded  the  threshold  value  for  high  gain.  In  many  cases  the  background 
noise  level  also  exceeded  the  medium-gain  threshold. 

37-  The  assumption  was  that  the  detection  ranges  predicted  with 
the  analog  computer  would  have  higher  values  and  exhibit  less  scatter 
than  those  measured  in  the  field.  The  measured  analog  records  were 
obtained  with  carefully  controlled  experimental  procedures  designed  to 
minimize  the  effects  of  variation  in  individual  footsteps  and  walking 
speeds  as  well  as  the  energy  loss  from  the  soil  to  the  geophone.  The 
detection  range  measured  with  the  DSID  more  closely  simulated  field 
deployment  conditions.  Therefore,  only  occasionally  would  the  DSID- 
measured  detection  range  be  as  great  as  those  predicted  with  the  analog 
computer.  The  detection  ranges  obtained  from  the  computer  were  plotted 
versus  measured  detection  distance  for  low  gain  (fig.  4).  In  most 
cases  the  computer-predicted  detection  ranges  compared  well  with  the 
maximum  measured  detection  ranges  and,  with  one  exception,  were  greater 
than  the  average  detection  range.  In  12  of  1 6  comparisons  the  com¬ 
puter  predictions  were  within.  5  m  of  the  maximum  of  the  measured  data. 

38.  Perhaps  more  significant  than  the  approximate  one-to-one 
correspondence  of  the  predicted  and  maximum  measured  data  is  the  fact 
that  the  range  of  distance  obtained  by  both  methods  correspond,  i.e. 

5-30  ra-  for  the  predicted  and  4-27  m  for  the  measured.  This  suggests 
that  both  procedures  are  about  equal  in  sensitivity  to  variations  in 
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COMPUTER  PREDICTED  DETECTION  DISTANCE,  METERS 

Fig.  4.  Predicted  analog  vs  measured  detection  ranges 

site  conditions.  On  the  "basis  of  this  analysis,  it  was  assumed  that 
the  analog  computer  predictions  could  he  used  to  estimate  the  quality 
of  predictions  obtained  from  the  digital  computer  program. 

39-  The  detection  ranges  determined  with  the  analog  computer 
were  used  as  a  reference  for  evaluating  the  digital  computer  program 
for  predicting  seismic  sensor  performance.  The  detection  range  values 
for  medium  and  low  gain  obtained  with  the  analog  computer  were  plotted 
against  the  low-  and  medium-gain  detection  ranges  predicted  by  the 
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digital  program  for  the  selected  test  sites  in  Panama  and  Arizona 
(figs.  5  and  6,  respectively). 

40.  The  curves  in  figs.  5  and  6  show  that  the  techniques  used  for 
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DIGITAL  PREDICTION,  METERS 


Fig.  5.  Comparison  of  analog  and  digital  predictions  of 
sensor  performance,  low  gain 


ANALOG  PREDICTION,  METERS 


predictions  cay  be  due  to  several  factors.  Firs*,  the  sensor  electronic 
components  look  at  the  entire  seismic  signal.  Thus,  the  peak  particle 
velocity  alone  does  not  convey  enough  inf  creation  regarding  the  charac¬ 
teristics  of  the  seismic  signal.  Second,  the  frequency  response  of 
the  sensor  and  the  frequency  characteristics  of  the  seissic  signal  have 
been  ignored,  and  these  factors  are  believed  to  be  very  significant  in 
predicting  sensor  performance. 
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PART  III.  TERRAIN  FACTOR  INPUTS 


1»1.  The  majority  of  the  terrain  factors  measured  during  the  field 
data  collection  were  analyzed  to  determine  the  highest  correlation  with 
the  peak  particle  velocity.  For  the  terrain  analyst,  however,  the  best 
model  may  be  of  little  value  if  the  input  requirements  are  so  stringent 
that  they  are  virtually  unattainable.  The  equation  used  for  the  digital 
computer  sensor  performance  model  was  selected  on  the  basis  of  the  rela¬ 
tive  ease  of  acquisition  of  the  input  values. 

1*2.  A  total  of  nine  terrain  factors  are  required  for  the  present 
version  of  the  model.  These  terrain  factors  are  listed  in  paragraph  26. 

The  ideal  way  to  obtain  the  input  values  is  by  direct  measurement  or 
calculation.  Often  access  to  the  ground  is  denied  and  the  terrain  fac¬ 
tor  inputs  will  have  to  be  estimated.  This  estimation  will  often  have 
to  be  made  from  aerial  photographs;  various  soils,  geologic,  physio¬ 
graphic,  or  land  use  maps;  and  other  literature.  This  part  of  the  re¬ 
port  presents  data  that  can  be  used  as  an  aid  in  estimating  the  terrain 
factor  value.  This  information  is  not  complete,  and  the  terrain  analyst 
should  supplement  it  with  information  from  other  sources. 

Estimating  Cone  Index,  Soil  Type, 

and  Water  Content  Values  J 

1*3.  The  cone  index  value  of  a  soil  is  a  function  of  soil  type 
(Unified  Soil  Classification  System  (USCS))  and  soil  moisture  content. 

The  type  of  the  surface  soil  can  often  b^  estimated  reliably  from  aerial 
photographs  or  soil  maps  (see  reference  6).  In  general,  the  water  con¬ 
tent  of  a  soil  cannot  be  readily  estimated  from  remote  sensing  imagery, 
visible  photography,  or  most  soils  maps  since  it  is  a  function  of  rain¬ 
fall,  topographic  position,  soil  type,  and  other  factors  related  in  a 
complex  manner.  Automation  of  a  soil  moisture  prediction  model  being 
developed  at  the  WES  is  scheduled  for  completion  in  June  1972,  thus  pro¬ 
viding  a  means  for  estimating  average  daily  soil  moisture  content.  If 
soil  type  and  water  content  are  known,  cone  index  values  for  the  0-  to 
15-cm  layer  can  be  estimated  from  the  generalized  relations  among  soil 


type,  moisture  content,  and  cone  index  shown  in  fig.  7-  Additional 
data  on  the  relation  of  cone  index,  soil  type,  and  water  content  are 
given  in  reference  8. 


CONE  INDEX 


Fig.  7.  Generalized  relations  among  soil  type  (after  reference  7), 
soil  moisture  content,  and  cone  index 

44.  Aids  for  estimating  surface  cone  index  are  not  readily  avail¬ 
able,  but  probably  could  be  developed  from  existing  data  at  the  WES. 
Surface  cone  index  values  are  often  30  to  50$  of  the  average  cone  index 
of  the  0-  to  15-cm  layer,  but  values  10$  of  the  cone  index  of  the  0-  to 
15-cm  layer  are  common. 
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Estimating  Dry  Density  and  Layer  Thickness  Values 


Mean  Dry 

uses 

Density 

Type 

g/cm3 

SM-SC 

1.62 

SC 

1.60 

SP-SM 

1.57 

SM 

1.50 

CL-ML 

1.49 

CL 

1.46 

ML 

1.37 

CH 

1.36 

OL 

1.32 

MH 

1.11 

OH 

1.00 

These  density  values  were  derived  from  analysis  of  approximately  1300 
samples  taken  from  the  15-  to  30-cm  layer  in  the  temperate  zone.  Nor¬ 
mally  the  0-  to  15-cm  layer  contains  more  organic  matter  and  will  ex¬ 
hibit  slightly  less  dry  density,  i.e.  values  from  10  to  15$  less  than 
those  shown.  Fig.  8  is  presented  to  allow  estimation  of  dry  density 

Q 

from  known  in  situ  moisture  content. 

Thickness  of  the  first  soil  layer 

46.  This  terrain  factor  value  is  probably  the  most  difficult  to 
obtain  by  noncontact  means.  Estimates  of  soil  thickness  can  often  be 
made  by  photo  interpretation  and  study  of  the  geologic  land  use ,  soil 
maps,  and  related  literature  of  the  area.  As  previously  mentioned, 
however,  these  estimates  may  not  be  valid  for  seismic  layers  that  are 
based  on  physical  soil  properties. 


Estimating  Comnression  Wave  Velocity 


47.  The  terrain  analyst,  having  determined  soil  or  rock  type,  can 
obtain  estimates  of  compression  wave  velocities  from  the  following 
tabulation . 
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GM/CM3  =  62.422  PCF 


MOISTURE  CONTENT  IN  PERCENT 

Fig.  8.  Dry  density  as  a  function  of  field  moisture  content,  0-  to  15-cm  layer 


Estimated  Compres¬ 
sion  Wave  Velocity 
_ m/sec _ 


Soil  or  Rock  Types 


Dry,  loose  topsoils  and  silts 

180  to  365 

Dry  sands,  loams,  and  slightly  sandy  or  gravelly 
soft  clays 

300  to  485 

Dry  gravels;  moist  sandy  and  gravelly  soils 

450  to  910 

Dry,  heavy  gravelly  clay;  moist  heavy  clays; 
cobbly  materials  with  considerable  sands  and 
fines;  soft  shales;  soft  or  weak  sandstones 

910  to  1450 

Water;  saturated  silts  or  clays;  wet  gravels 

1460  to  1524 

Compacted  moist  clays;  saturated  sands  and 
gravels;  soils  below  the  water  table;  dry  medium 
shales;  moderately  soft  sandstones;  weathered 
moist  shales  and  schists 

1460  to  1829 

Hardpan;  cemented  gravels;  hard  clay;  boulder 
till;  compact  cobbly  and  bouldery  materials; 
medium  to  moderately  hard  shales  and  sand¬ 
stones;  partially  decomposed  granites;  jointed 
and  fractured  hard  rocks 

1676  to  2b 3d 

Hard  shales  and  sandstones;  interbedded  shales 
and  sandstones;  slightly  fractured  limestones 
and  crystalline  rocks 

2438  to  3657 

Unweathered  limestones,  granites,  gneiss,  and 
other  dense  rocks 

3657  to  6100 

Rayleigh  Wave  Velocity 

48.  Rayleigh  wave  velocity  is  equivalent  to  shear  wave  velocity. 
For  very  rough  estimates,  Rayleigh  wave  velocity  can  be  assumed  to  be 
40  percent  of  the  compression  wave  velocity. 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


49.  Results  to  date  in  the  IVES  seismic  sensor  programs  indicate 
that  terrain  factors  can  be  correlated  with  the  seismic  response  charac¬ 
teristics  of  an  area.  These  correlations  can  then  be  used  to  predict 
the  peak  amplitude  of  a  seismic  signal  with  distance;  however,  predic¬ 
tion  of  seismic  sensor  performance  by  this  method  is  by  no  means  ade¬ 
quate.  Although  the  present  version  of  the  computer  program  is  inade¬ 
quate  for  predicting  SID  performance,  it  can  be  easily  updated,  and  an 
improved  prediction  capability  can  be  expected  as  additional  information 
becomes  available  and  additional  analyses  are  completed. 

50.  Although  the  predicted  peak  particle  velocity-distance  curves 
are  not  adequate  to  allow  accurate  determinations  of  sensor  performance, 
they  may  be  used  with  confidence  (within  the  range  of  experimental  data 
upon  which  they  are  based)  to  obtain  a  relative  comparison  of  the  seis¬ 
mic  response  characteristics  of  difference  areas. 

51.  The  terrain  factor  values  required  as  inputs  to  the  computer¬ 
ized  model  are  those  that  are  common  to  the  earth  sciences.  Consider¬ 
able  data  in  various  parts  of  the  world  have  been  collected  on  these 
terrain  factors,  and  reasonable  estimates  of  their  values  can  be  made 
on  the  basis  'Of  literature  and  other  information  sources,  such  as  re¬ 
mote  sensing  products  and  the  WES  soil  moisture  prediction  system. 

Recommendations 


52.  Continued  work  on  the  development  of  an  analytical  procedure 
for  predicting  SID  performance  is  recommended.  Emphasis  should  be 
placed  on  developing  theoretical  equations  for  predicting  the  complete 
particle  velocity  wave  train  as  a  function  of  distance.  Computer 
routines  capable  of  .interpreting  the  particle  velocity  wave  train  in 
a  manner  analogous  to  that  of  the  sensor  should  be  -developed.  These 
sensor  simulation  routines  should  be  combined  with  the  theoretical  wave 
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propagation  equations  to  provide  a  theoretically  based  SID  performance 
model . 

53-  Additional  research  should  be  conducted  to  develop  improved 
techniques  for  obtaining  the  terrain  factor  input  values. 
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APPENDIX  A:  COMPUTER  PROGRAM  TO  PREDICT  SID  DETECTION  DISTANCE*' 


Description  of  Function  and  Method 

1.  A  computer  program  has  been  developed  to  estimate  detection  dis¬ 
tance  (in  the  program  detection  distance  is  termed  detection  radii)  from 
one  man  walking  to  a  seismic  sensor  as  a  function  of  soil  parameters .  The 
program  is  intended  to  be  self-instructive,  easy  to  learn  to  use,  and,  since 
it  operates  in  the  conversational  mode,  suited  to  processing  of  relatively 
small  quantities  of  data  input  from  a  teletype.  The  peak  particle  velocity- 
distance  equation,  described  in  the  main  text,  is  computed  in  a  subroutine 
so  that  it  can  be  revised  with  little  change  in  the  main  program. 

2.  The  values  of  the  peak  particle  velocity  (resulting  from  a  hammer 
drop*-*)  are  computed  for  radii  of  2,  5,  10,  15,  20,  30,  and  ko  m.  If  the 
computed  peak  particle  velocity  becomes  zero  or  negative  at  a  radius  of 
less  than  40  m,  velocity  values  corresponding  with  greater  radii  are  ig¬ 
nored.  The  array  of  values  for  radius  vs  peak  particle  velocity  is  entered 
into  a  spline  curve-fitting  routine  (subroutine  SPL),  which  generates  a  peak 
particle  velocity  for  each  meter  of  radius  from  2  to  40  m.  These  values  are 
written  into  a  disc  file  (fig.  Al). 

3.  After  the  spline  routine  is  completed,  the  resulting  array  enters 
a  subroutine  called  PRED,  in  which  the  peak  particle  velocities  required  to 
open  the  sensor  logic  at  medium-  and  low-gain  settings  are  found.  The  ra¬ 
dius  associated  with  the  selected  peak  particle  velocity  becomes  the  pre¬ 
dicted  detection  radius.  If  the  associated  radius  exceeds  40  m,  the  solu¬ 
tion  is  not  valid,  and  the  value  of  the  associated  radius  is  set  to  zero. 


*  This  program  is  furnished  by  the  Government  and  is  accepted  and  used  by 
the  recipient  with  the  express  understanding  that  the  United  States 
Government  makes  no  warranties,  expressed  or  implied,  concerning  the  ac¬ 
curacy,  completeness,  reliability,  usability,  or  suitability  for  any 
particular  purpose  of  the  information  and  data  contained  in  this  program 
or  furnished  in  connection  therewith,  and  the  United  States  shall  be 
under  no  liability  whatsoever  to  any  person  by  reason  of  any  use  made 
thereof.  The  program  belongs  to  the  Government.  Therefore,  the  re¬ 
cipient  further  agrees  not  to  assert  any  proprietary  rights  therein  or 
to  represent  this  program  to  anyone  as  other  than  a  Government  program. 

**  Equivalent  to  a  footstep. 


Al 


Inrnt  variables 


Vari¬ 

able 


Record 

Name 

Column(  s ) 

,  Mode 

Units 

Description 

1 

C115 

Free-fi°ld 

Floating 

Dimensionless 

Cone  index 

(0-  to  15- cm  layer) 

wc 

Free-field 

Floating 

Dimensionless 

Water  content  of 
surface 

WC1 

Free-field 

Floating 

Dimensionless 

Water  content  of 
first  layer 

VP 

Free-field 

Floating 

m/sec 

Compression  wave 
velocity 

VR 

Free-field 

Floating 

m/sec 

Rayleigh  wave 
velocity 

H 

Free-field 

Floating 

m  . 

Thickness  of  first 
seismic  layer 

EFS 

Free-field 

Floating 

Dimensionless 

Percent  fines 
decimal 

GAMD 

Free-field 

Floating 

gm/cc 

Dry  density  of 
surface 

GAMD1 

Free-field 

Floating 

gm/cc 

Dry  density  of 
first  layer 

2 

ANAM 

1-6 

A 

Site  identification 
alphanumeric 
characters 

3 

NAME 

1-6 

A 

-- 

Name  of  saved 

file  into  which 
to  write  results 


4 

XO 

Free-field 

Floating 

m 

of  equation 

Initial  radius 

XMAX 

Free-field 

Floating 

m 

Maximum  radius 

DELX 

Free-field 

Floating 

m 

Spline  fit  increment 

• 

desired 

Output 

4 .  An  example  of  the  output  from  program  KLN012  with  teletype 
input  of  soil  parameters  is  shown  in  fig.  A2.  Circled,  numbers  in  the 
left  margin  correspond  to  the  input  record  order  in  the  input  description 
above . 


5.  Save  the  names  of  the  output  files.  If  no  plots  of  radius 


I 


f  / 


versus  peak  particle  velocity  are  desired,  the  same  name  may  "be  used  for 
several  sites. 

6.  Run  the  program.  Data  must  be  entered  from  the  keyboard  as 
shown  in  fig.  A2  when  it  is  requested  by  the  program. 

Program  Listings 


7.  The  program  listing  is  shown  in  fig.  A3.  The  language  used  is 
FORTRAN  IV  adapted  for  use  on  a  G-^37  (Honeywell)  time-sharing  computer. 

Flow  Charts 

8.  The  flow  charts  of  the  main  program  and  necessary  subroutines 
are  presented  in  plates  A1  and  A2. 


A3 


4 


} 
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NAM1 

DISTANCE 

PARTICLE 

VELOCITY 

INCREMENT  (M) 

(CM/SEC  x  1(T3) 

00001 

TEST' 

00002 

2.000 

■  12.743 

00003 

3.000 

10.438 

00004 

4.000 

8.307 

00005 

5.000 

6.522 

00005 

6.000 

5 . 268 

00007 

7.000 

4.424 

000  OS 

'  3.000 

3.830 

000  OS 

9.000 

3 . 5.26 

00010 

•  10.000 

3.253 

0001  1 

1 1.000 

2.979' 

00012 

12.000 

2.690 

00013 

13.000 

2-.  403 

000  M 

14.000 

2.130 

00015 

15.000 

•  1.389 

0001 5 

16.000 

1.687 

COO  17 

17.000 

1.523 

000  is 

18.000 

.  1.389  ■ 

000  is 

19.000 

1 . 27S 

00  0  20 

20.000 

1.173; 

00021 

"  '  21.000 

.  1.037 

00  022 

22.000 

1.003 

00023 

23.000 

0.925 

00024 

24.000 

0.352 

00025 

25.000 

0.785 

000  25 

26.000 

0.724 

00027 

27.000 

0.668 

00028 

23.000 

0.617 

00.029 

29.000 

0,570 

00030 

o  0  *  0  0  0 

0.529 

00031 

31,000 

0.491 

00032 

'  32.000 

0.458 

00033 

33;000 

0.429 

00034  ' 

34,000 

0.402 

00035 

35.000 

0.378 

00035 

36,000 

0.357 

00037 

37,000 

0.337 

00038  ' 

3S.000 

0.319 

00039 

39,.  000 

0.301 

00040 

40,'.  000 

0,284 

Fig.  Al.  Peak  particle  velocity  values  generated  by  the 
spline  curve-fitting  routine 


RUN 


XLN012  14:23  WES  01/2G/72 


ZNTEK  STOP 'IF  YOU  WISH  TO  STOP  THE  PROGRAM 

INPUT  OS  FIRST  LINE,  SEPARATING  THE  VARIABLES  BY  COMAS , 

THE  FOLLOWING  VARIABLES.: ( 1  )CONE  INDEX (6- 15  CW), 

(2) WAT EH  CONTENT  Of  SURFACE,  (3 )WAT£R -CONTENT  OF  FIRST  LAYER  , 
(^COMPRESSION  WAVE  VELOCITY  (ir/SZC ) ,  (5  )RAYLEI  GH  WAVE  VELOCITY  CM /SEC )  , 
(6)  THICK  NESS  OF  FIRST  SEISMIC  LAYER  (CM) ,  (7)  PERCENT  FINES  (7.), 

(8)  DRY  DENSITY  OF  SURFACE  (CM /CC  )  ,  (9)  DRY  DENSITY  OF  FIRST  LAYER  (Gil /CC ) 

.  ON  SECOND  LINE  ENTER  UP  TO  SIX  CHARACTERS  3F  IDENTIFICATION 


0  7120.,. 38.. 45; 35 1 . , 134. ,488. ,82. ,.82,lil 

(D  ?test  ' 

^ENTER .FILE  NAME  INTO  WHICH  TC  WRITE  SPLINE  VALUES 
Qjy  ?NAi'i  1 

COMPUTED  VALUES  OF  R  VS.  UDH 
SITE  TEST 


R 

UDH 

2.000 

•  12.743 

1 

5.000 

'  6.522 

2 

10.000 

3.253 

3- 

15.000 

1  .889 

4 

20.000 

I".  178 

5 

.30.000 

0.529 

6 

40.000 

0.284  . 

.  7 

ENTER  INITIAL  R  AND  MAXIMUM  R  FROM  ABOVE  TABLE  .AND 
INCREMENT  OF  SPLINE  VALUES  DESIRED(USUALLY  1) 

@  72,40,1 

DETECTION  RADII  FOR  SITE  TE 
MEDIUM  LOW  CAIN' 

0.  37. 


Fig.  A2.  Sample  run  of  program  KLN012 


KLfJC  1 2 


r-T  ry ,  1  20 

J!:'r  T;:iS,  P*£C[lAK  US£S  c‘^  G"  TWO  ECUAlIOr.S  TO  PREDICT  JETECTIGk  RADII 
120C  Cr  A  Sf.r,SCR  AT  THREE  DIFFERENT  GAIN  SETTINGS  ■ 

13G  DIi'.ERSI Oi.  ,((7)  .UDH  (7)  ,i'DH2(7)  ■ 

135  EQUI  VALENCE  (UDH, UDH  2) 

140  tiz.nL  li A i‘J w 

15  OC  THESE  DATA  ARE  THE  RADEI  FOR  WHICH  UDH  VALUES  ARE  COMPUTED 
160  DATA  R/2., 5., 10. ,15. ,20. ,30. ,40./ 

165  DATA  TEMP/"  "/  ' 

130  10  COS.t!  HUE 

240  51  FORMAT (  INPUT  Or*  FIRST  LINE,  SEPARATING  THE  VARIABLES  BY  CO/i.-iAS , 
25 03./ "THE  FOLLOW! KG  VARIABLES :<1  )CCKE  INDZX(0-15  CM), 

252*"/“ 

255A(2HA7Et?  CONTENT  CK  SURFACE,  (3JWATER  CONTENT 
SCODOF  FIRST  LAYER  ,  "/“( ^COMPRESSION  W&U?  VELOCITY  CM/SEC ) . 

2653  (5 )RAYLEI GH  WAVE  VELOCITY  (K/SEC ) , "/"(6>THICKNESS  OF  FIRST  SEISMIC 
2553 

2703LAYER  (CM  ) ,  (  7)  PERCENT  FI  RES  (7. ) ,  "/  (3)  DRY  DENSITY  OF  ,, 

2303  SUKFACETGM /CC  )  , (9)  DRY  DENSITY  OF  FIRST  LAYER (fi;-i/CC) 

2353//"0.’-  SECOND  LIRE  ENTER  UP  TC  SIX  CHARACTERS  OF  IDENTIFICATION  //) 
270  PRINT, “ENTER  STOP  IF  Y.CU  WISH  TO  STOP  THE  PROGRAM 

§5'  295  PRINT  51  ’ 

30 0  149  READ  ,CI  15  ,WC  ,VC l ,VP ,VR ,H ,PFS ,GAKD ,GAMD1 
305  READ  3 06, A HAM 
5C5  306  FORMAT TAG) 

512  I  P(At.Art  ,EC  .  6HSTOP  )CALL  EXIT 

320  320  PRINT , "ENTER  FILE  NAME  INTO  WHICH  TO  WRITE  SPLINE  VALUES" 

350  READ  1 5 4,  NAME 
340  1 54  FORMAT (A 5) 

400  30  CALL  THKE£(CI  15  ,WC  ,VC'I  ,VP  ,VR ,H  ,PFS ,GAMD,GAMD l  ,UDH  ,R ) 

42 OC  SOLVE  FOR  SELECTED  VALUES  CF  UDH 
550  X1=J-1 

570  160  DC  250  1=1,7 
530  I  F  (UDH 2(1  ) )  l~70 ,170,250 
590  170  K2=T 
600  UDH2(I )=0. 

51  OX  267  FORMAT ( 2(1 5, 2F1 0.2)) 

620  GOTO  250 
630  250  CONTINUE 
640  K 2=1-1 
545  260  CONTINUE 
660  660  COKTPNUE 

570  203  /FORMAT ( "COMPUTED  VALUES  OF  R  VS.  UDH  7 
6503 2X  i "SI TE ",2X. A6//.7X  ,  "R " , 6X  ,  "U Dil "/ ( 2F  10.3,15)) 

63  0  PR  I N  T  233 ,  AKAH  ,  (R  (IZ  )  ,UDH  (I  Z ) ,  IZ  ,  I  Z=1  ‘.K  2) 

720  293  CALL  SPL(K2,R  ,UDH  ,3  NAl’iE  ,ANAH  ,KNT2) 

730  CALL  PRED(3,NA»'.£,KkT2)  '  ' 

740  GO  TO  10 


Fig.  A3.  listing  of  program  KLN012 
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XLNC‘12  CONTINUED 


750  r-S3  STOP  3 S3 
750  END  ■ 

1030  SUBROUTINE  RRED(NDEV  .NAME  ,Ki\T) 

1030C  THIS  RC-UFli.’E  CGE-PAKES  VALUES  FOR  LCa'j 
U00C  TO  VALUES  OK  THE  R  VS.  UDH  CURVE. 
1110  DI.'iEivSICii  R ( 100), UD( 100), FRC3) 

1120  DIMENSION  DRAD(3)  . 

1130  DATA  FR/.l,  .2,0.5/ 

1140  CALL  OPENF  ( NDSV  .NAME ) 

1150  READCNDEV ,825)  ,ANAM 
1160  S26  FORMAT  (A  6 
1 1  705) 

1130  DO  662  1=2,3 
HOC  562  DEAD ( I )=0. 

1200  READ  (N’DEV  ,329)  (R  (I  )  ,UD(I  )  ,1  =  1  ,K11T) 

1210  829  FORMAT (2F1 0.3) 

1220  DC  ! 50  '  =2  3 

1230  DC  50J=  1  ,K!)T 

1240  I F(FE(L)-UD<J))50, 100,100 

125 C  50  CONTINUE 


m£DI  lii'i , 


AiO  HIwjH  GAIk 


1270  GO  TO  150 
1230  100  DRAD (L  )=R  <J ) 

1290  150  CONTI  RUE 

1295  PRINT  1296, AMAH 

1296  1296  FORMAT  ( 10X  ,  "DETECTION  RADII  FOR  SITE”,  IX  ,A2) 

1300  PRINT  1301  . 

1301  1301  FORiiAT (13X,  "MEDIUM ",i OX, "LOW  GAIN") 

1310  PRINT  131 1 , ( DR A 0 (K ) , K  =  2 , 3 ) 

1311  1311  FCRi*.A,T-(5X,2F13.0) 

1320  PRINT  132-1 

1321  1321  FORMAT (////) 

1330  CALL  CLCSEr(NDEV) 

1340  790  RETURN 

1350  END  . 

1350  SUBROUTINE  SPL  (B  ,X  ,Y  ,NDEV  jNAME  ,ANAM  ,XNT ) 

1370C  SPLINE  TEST  PROGRAM  BY  J  CHEEK,  ADPC 
1330  DIMENSION  X ( 100 ).,Y  <  1 00) , TORQUE ( 1,00) 

1590XPSI  NT  836  ,NDEV  ,KA,‘.E 
1400  336  FORMAT (1 5 , AS) 

1410  CALL  0PE1IF  (.“DEV, NAME) 

1420  KuT=0 

1430  10  CONTINUE 

1440  CALL  SPLINE  (X,  Y,  K,  TORQUE) 

1450XPRI NT , "WE  ARE  NOV  IN  THE  SPLINE  ROUTINE”  # 

1460  PRINT, "ENTER  INITIAL  R  AND  MAXIMUM  R  FROM  ABOVE  TABLE, AND 
1465  PRINT ."INCREMENT  CF  SPLINE  VALUES  DESI RED (USUALLY  1) 

1470-  WRITE^DEV.,  157), ANAM.  ‘ 

1430  157  FORMAT (A 6) 

■1490  READ,  XO.XMAX,  DELX 

1500  100  CONTINUE 
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•f  i 
t  ■ 

ft 


OLY  -  YL)  /  H  2ZZZ 


KLKC  1  2  CO.’lTlSfiED 


1530  KKT=KM+1 

154.C  120  FC3.-,AT(2F10.3)  . 

1550  X0  :  Xil  +  9ELX 
1560  IF (XO.GT.Xi'iAX)  GO  TO  110 
1570  GO  TO  1 00 
1530  110  CONTINUE 
1550  CALL  CLOSER (rJDZV) 

1600  /tETURN 
16 1 0  EN3 

1520  SUBnOUTI  NE5PLTRE  (XI,Y , N  ,S2) 

1630  D1  MEllSlOli  ,XCl)  ,Y(  1),S2(-I ) 

1640  DATASPSLB/-1./  ' 

1650  K1  s  s  -  l 
1560  ASSI GN54T0ISU 
.1670  DC59I  =  1  ,f;j  ‘ 

1530  H=xa+i>-xci ) 

1630  0LY=  CY  (I  +1 (I  ))/H 
1700  GO  101  b>  ,  04  ,  33; 

1710  53i(2Z-ZZ=l<L +K 
1720  52(10  =  2.  * 

1730  GCT055 
1740  54ASSiGN'53TOISW 
1750  551(L=1!  '  . 

1760  YL=3LY 
1.770  59CCBT1BUE 
.173  OC 

1790  S2(<1  )  =  0, 

1800  S?.(!»=0. 

13,10  OMEGA  =  -1.0.7 17968 
1820  5Z  TA  :  0 , 

1330  ASSIGN  154  TO 
1840  00  10  I  =  1  N 1 

1350  H  =X'(  I  +1~)  -X  CI  ) 

1850  DLY-CY  CI+1,)-Y  (I ) )  /H 
IS 70  GO  TO.  ISWi ,( 154,  153) 

1330  153H2ZZZ=HL;+H- 
1690  81=  .  5*HL./H2ZZZ 

1900  y  =  Col  *  $2(1-1.)  +  (.5  -  Cl)  *  S  2(1  +  1 )  +  S2(I) 

;  Sm  1 *  oiCGa 

1930  Z  =  4CS (W) 

1940  154- ASSIGN  153  70  l Sill 
1950  156  ilL=M 
1950  YL=9LY 

1970 '10  CONTINUE  , ' 

193  0  I  f'(  ETA  -  EfiSLii  914,5,5 
1990  1;4  C01JTI rfl)E 
2000  HETU.TN 


ISWI 


I 

1 


!  i 


i 


jjl 


< 

1 
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KLtJO  1 2  COI.TI  SUED 


1 0  =.ND' 


2020  Sur'SCUTIKESFLlM  (XX  ,FXX  ,rfXX  .FPrXX  , 


2u30.’<X  ,Y,5-i,i;  ,,‘i  ,16) 

20-50  oi:-i£;:sic:.x(i),Y(i)ts2(i) 

2050  0 

20 SO  XP=XX 

20  70  1  =  1 

2030  IF(XF-X(i) >52,17,55 
209  0  5  2.v;=- 1 
2100  XP=X(1) 

2110  ROTO  17 

21 20  551 F (XP-X (K ) )57 ,59,53 
2130  5SIr(XP-X(I  ))60, 17,57 

21  AO  571  =  1+1 
2150  RO  TO  5S 

21  so  58  1. 

2170  XP=X(R> 

21 5  0  39  I  -  ii 
2190  601=1-1 
2200  1 -7H  T  1=XF-X(I ) 

£210  HT2=XP-X(I+1) 

2220  PROD:. IT  1*1112 
2230  DX=X(I+1)-X(I ) 

2240  DELY=(Y(I+1)-Y(I  ))/DX 
2250  S  3=  (S  2(1+1  )-S2(I  ))/DX 
‘226  0  FPPXX=5  2(I  )+HTl*S3 

22  7  0  DELSCS :  (S  2(1  )+S2(I +D+PPPXX  )/Sy 
2230  FXX=Y (I >+H7 1 *DELY  +PROD*DELSQS 

2290  FPXX  =  DELY  +  (HT 1  +  H T 2):  *  DELSGS  +  PROD  *  S3  /  S. 

2500  I F (F ,£Q, 0)R0TC 1 00 
2310  FXX  =  FXX-I-FPXX*('XX-XP) 

23  20  1  00CC0TIMJE 
2330  KtTL’RN 

234  0  END 

4030  SULROUTIfiE  THREE  (Cl  l-5,,v.'C  ,V!C  l  ,VP  ,VR  ,il  ,P,FS  ,GAilD,GAi'lDl  ,UDH  ,R ) 
4005  DIBEKSIOK  UDH(l),R(t) 

40 1  0  A  1  =-S,02*GAi'iD+1  7. G5 *!1<*K / ( VP *VP ) 

4020  A  2=-  1 . 1  H5*il*ll*H/(  VP*VP*CI  1 5  ) 

4030  A  3 = - 4 # 7 6 * ( R A rt D *‘*’C / (4  .-GAf.D/2.G5)>*H**4/VP**4 
404  0  A  4  =  0 . 0 1 25/  ( 1 .  - VH  /VP )  *GAr; D *WC 7 ( 1 .  -CAHD/2. 65 ) 

4050  A  =  A  1 +A 2+ A  5+A  4+8,31 

4C50  A  1  =  . 000157*11 

4370  A  2=  .326* ( 1  »+WC  )*GAKD 

40-30  A  3= .  09 2*U  .  +WC  1 )  *GAK D 1 

4090  A4  =  -4.C-3*H*VR/PFS*WC  1 

4100  A5  =  -,00S2*GAl-lD*UC/(  1  ,-GA.'.D/2.  S5  > 

41 1 CX  *11  **4 /VP  **4 

41 20  A  1  =  A  1 +A 2+A 3+A 4+A5 - . 1 69 

4125  DO  4135  1=1,7 

4130  UDIKI  )  = ,  1 !+A*EXP  C-M*R  (I )  3/SSRT  (R(I  >> 
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KLK012  CONTINUED 


41 35  4135  GOiiTlilUE 
4193  RETURN 
4200  El.'D  • 
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SUBROUTINE  SPLINT 
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HT2)  ♦  OELSQS  ♦  PROD 
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111.  ABSTRACT 
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Improved  guidance  manuals  for  planning  the  deployment  and  emplacement  of  seismic  intrusion  detectors  (SID's) 
are  needed  to  optimize  the  use  of  these  devices  for  battlefield  surveillance.  The  development  of  thes-> 
Military  Geographic  Intelligence  (MGl)  products  requires  a  detailed  understanding  of  the  operating  princi¬ 
ples  of  the  detector  coupled  with  an  equally  detailed  understanding  of  the  interactions  of  the  sensor  prop¬ 
agation  mode  with  the  operational  environment.  This  report,  presents  the  results  of  a  preliminary  analysis 
of  data  collected  in  a  wide  range  of  environments  at  22  sites  in  Panama,  10  sites  in  Puerto  Rico,  6  sites 
near  Yuma  Proving  Ground,  Arizona,  and  9  sites  near  Ft.  Huachuca,  Arizona.  Multiple  regression  techniques 
were  used  to  determine  the  terrain  factors  that  could  be  correlated  with  the  seismic  responses  resulting 
from  a  man  walking  or  a  controlled  source  (drop  hammer)  that  simulated  the  signature  resulting  from  a  foot¬ 
step.  The  measure  of  seismic  response  was  peak  particle  velocity  as  a  function  of  distance  from  the  source. 
The  terrain  factors  that  correlated  best  with  peak  particle  velocity  were  the  thickness  of  the  first  refrac¬ 
tion  layer,  cone  index  of  the  0-  to  15-cm  soil  layer,  dry  density  of  surface  soil  and  first  soil  layer, 
water  content  of  surface  soil  and  first  soil  layer,  compression  wave  velocity,  Rayleigh  wave  velocity,  and 
grain-size  distribution.  An  empirical  equation  was  developed  to  predict  peak  particle  velocity  versus  dis¬ 
tance  as  a  function  of  the  terrain  factors.  The  particle  velocities  required  to  trigger  the  logic  of  the 
Phase  III  SID's  were  superimposed  on  the  predicted  peak  particle  velocity  curves  to  arrive  at  a  prediction 
of  sensor  performance.  These  computation  procedures  were  computerized  to  make  a  prediction  model  for  rela¬ 
tive  SID. performance  as  a  function  of  terrain  factor  values.  The  empirical  prediction  equation  adequately 
predicted  the  peak  partible  distance  relation;  however,  the  predictions  of  sensor  performance  were  inade¬ 
quate;  The  errors  in  the  predictions  of  sensor  performance  were  attributed  to  the  inadequacy  of  the  peak 
particle  velocity-distance  relation  to  represent  the  complex  interaction  of  the  entire  seismic  signal  and 
the  sensor.  Frequency, .characteristics  of  the  seismic  signal  and  the  frequency  response  characteristics  of 
the  sensors  also  must  be  considered. 
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